REACTION-DIFFUSION FRONT SPEED ENHANCEMENT BY FLOWS 



ANDREJ ZLATOS 

Abstract. We study flow-induced enhancement of the speed of pulsating travehng fronts 
for reaction-diffusion equations, and quenching of reaction by fluid flows. We prove, for 
periodic flows in two dimensions and any combustion-type reaction, that the front speed is 
proportional to the square root of the (homogenized) effective diffusivity of the flow. We 
show that this result does not hold in three and more dimensions. We also prove conjectures 
from [U [31 [H] for cellular flows, concerning the rate of speed-up of fronts and the minimal 
flow amplitude necessary to quench solutions with initial data of a flxed (large) size. 



1. Introduction and the Main Results 

It is well known that the presence of a fluid flow can significantly increase mixing properties 
of diffusion. The study of this phenomenon, sometimes called eddy diffusivity, has been 
the aim of a large body of mathematical and physical literature. Questions of long time- 
large scale behavior are usually addressed via techniques of homogenization theory (see, e.g., 
[l3l[2T] and references therein). This approach is appropriate when one can wait a long time 
for mixing to take effect. The presence of other processes, however, may introduce additional 
time scales to the problem. One such process is reactive combustion, which happens on short 
time scales and therefore requires a different approach to the study of combustive mixing. 

The effects of flows on combustion have recently been studied by various authors, both 
qualitatively and quantitatively [IiaEllIllHliaillllianSllISlinillHlESlESlETllSglH 
The main effects are two-fold. A strong flow can speed up propagation of a reaction (such as 
a wind spreading a fire) but also extinguish it (the "try to light a match in the wind" effect). 
The models used are reaction-advection-diffusion equations, in which the first phenomenon 
is manifested by the enhancement of speed of their (pulsating) traveling front solutions, and 
the second by quenching of solutions with (large) compactly supported initial data. 

In the present paper we study both these effects for stationary periodic flows. We consider 
the reaction-advection-diffusion equation 

Tt + u{x)-VT = AT + f{T) (1.1) 

for the (normalized) temperature T{t,x) G [0, 1] of a premixed combustible gas, with {t,x) G 
M X M.'^. The gas is advected by a periodic incompressible (i.e., V ■ u = 0) mean- zero vector 
field u G C^''^(]R^) (also calld flow). For the sake of simplicity, we will assume that all periods 
of u are 1, that is, the cell of periodicity is the d-dimensional torus T'^ = [— |, |]'^, with — | 
and I identified. The general periodic case is identical. The non-negative reaction function 
f G C^'^([0, 1]) accounts for the increase of temperature due to a chemical reaction such as 
burning, and is of the combustion type. That is, there is ^ G [0, 1) such that f{s) = for 
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s G [0, 9] U {1}, /(s) > for s G {9, 1), and / is non-increasing on (1 — 6, 1) for some 6 > 0. If 
9 > 0, then / is an ignition reaction (with ignition temperature 9), otherwise / is a positive 
reaction. A special case of the latter is the Kolmogorov-Petrovskii-Piskunov (KPP) reaction 
[I9j with < f{s) < sf (0) for all s G (0, 1). 

A pulsating traveling front in the direction of a unit vector e G M*^ is a solution of (11.11) of 
the form T(t, x) = U{x ■ e — ct, x), with c G M the front speed and [/ : M x R'^ — > [0,1] periodic 
in the second variable such that 

lim U{s,x) = 1 and lim U{s,x) = 0, (1.2) 

S—f — OO s— > + oo 

uniformly in x G M'^. It is well known that under our assumptions on u and /, for each e G 
there is a unique c*{u, /) > such that a pulsating traveling front in direction e and with 
speed c exists if and only if c = c*{u, f) for ignition reactions [28], resp. c G [c*(m, /), oo) for 
positive reactions [2]. 

In both cases we will be interested in the fronts with the unique /minimal speeds cl{u, /). 
These are the most physical ones because they determine the speed of spreading of solutions 
to the Cauchy problem for (II. ip with (large enough) compactly supported initial data. An 
exact formula for €*(«, /) has been obtained for general reactions in [TU] (and for the special 
case of KPP reactions also earlier in [1]). Unfortunately, it is a complicated variational 
expression and it is not obvious how to use it to obtain simple general estimates on c*{u, /). 

Our goal is to derive simple estimates on the front speed c*(m, /), and use them to answer 
open questions from [H [3l |TT] concerning speed-up of pulsating fronts and quenching of 
reaction by strong flows. We will provide such estimates in two dimensions, in terms of the 
size of / and the effective diffusivity Df.{u) of the flow u (in the direction e). The latter 
quantity can be found in a much simpler way than c*(u,/) using (II. 5p and (11.60 below. It 
appears in the homogenization theory which, as mentioned above, is applicable to the study 
of long time behavior of the solutions of the related linear PDE 

i)t + u{x) -S/i) = Ai). (1.3) 

Despite the fact that the presence of reaction introduces a short time scale to the model, 
we will be able to show by other methods that in two dimensions (but not in three!), the 
effective diffusivity determines the front speed up to a bounded factor. 

The long time behavior of the solutions of (II. 3p is governed by the effective diffusion 
equation 

^t = V-{D{u)V^), (1.4) 

where the (x- independent positive symmetric) effective diffusivity matrix D{u) is obtained 
as follows. For any e G M'^, let Xe{x) be the periodic mean-zero solution of the cell problem 

- Axe + M ■ Vxe = ■ e (1.5) 

on T'^. Then D{u) is given by 



e • D{u)e = / (Vxe + e) ■ (Vxe' + e')dx = e - e + / Vxe ■ VXe'dx 
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for any e, e' G M*^. The effective spreading in the direction of a unit vector e G M"^ is now 
governed by the effective diffusivity 

D,{u) = e ■ D{u)e = 1 + || Vxe|li2(T.). (1.6) 

When the nonhnearity in (11.11) is weak (the reaction time scale is large) so that we have 

Tt + u{x)-VT = + ef{T) 

with £ <^ 1, the long time-large space scaling t^t/e^.x^x/e gives 

Tt + \u -VT = AT + f{T). 

The homogenized version of this equation is 

Tt = V-{D{u)VT) + f{T), 

with the unique/minimal front speed in direction e depending on De{u) and /. Although 
the above approximation holds only on certain space-time scales, it does suggest a relation 
between c*(m, /) on one hand and Df,{u) and / on the other. Our main result confirms this 
relation in two dimensions: 

Theorem 1.1. Let f he any combustion-type reaction. There are Ci(/),C2(/) > such that 
for any 1-periodic incompressible mean-zero C^'^ flow u on and any unit vector e G M^, 



C,{f)^DM<c:{uJ) < C2{f)^DM- (1-7) 

Remarks. 1. We have the following estimates on the constants in (II. 7p . From the results 
of [25] and monotonicity of c*(u, /) in / it follows that 



C2U) < C^\\f{s)/s\\^ ( 1 + ^\\f{s)/s 



for some C > 0. If m(^{f) = min{/(s) | s G [C, 1 — ]} > (for each / as above there is 
such ( G (0, 1)), then from the proof of Theorem II. II it follows that 

for some > 0. We note that this estimate is optimal up to a constant for small m(^{f) 
(due to (II. Sp ). but also for large mc_{f). Indeed, Ci(/) is uniformly bounded above in / for 
KPP reactions |25j and thus for all reactions. 

2. In particular, Cj{mf) ~ {j = 1, 2) as m — > for any fixed /. 

3. In the case of KPP reactions, this result has been proved in [231 125]. This case is 
considerably simpler due to the fact that c*(m, /) = c*(m, /), where /(s) = /'(0)s and c*(m, /) 
is the minimal front speed for the linear equation (ll.ip with / in place of /. (Fronts for / do 
not converge to 1 as a; ■ e ^ — 00 but rather grow exponentially.) 

The relation cl{u, f) ~ C{f)^e (u) is analogous to that in the case of m = and constant 
diffusivity matrix D > 0. Indeed, spatial scaling x t— > \/Dx shows that the unique/minimal 
front speed in the direction e for Tt = V-{DVT) + f{T) is C{f)Ve-De (with C{f) = c*{0, f) 
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independent of e). The problem with u ^ 0, however, is that the convergence of solutions of 
f ll.3l) to those of fll.4l) occurs on large time scales, while solutions of fll.ip can be effectively 
estimated using (11.31) on short time scales only. We will therefore study short time diffusivity 
for (11.31) in Section [2] and relate it to De{u) in Theorem 12.11 below. The theorem, which 
applies in all dimensions, will be a key step towards overcoming this problem for d = 2. 

The restriction to two dimensions comes from Theorem 13.11 below, a relation between the 
speed and the width of a front. It turns out, in fact, that this result and Theorem 11.11 are 
false in dimensions d>3, and Theorem 15.11 shows that the bounds in (I1.7P cannot hold with 
flow-independent Ci(/) and C2(/) for d> 3\ 

This raises an interesting question about large deviations of the stochastic process Xf 
from (12.11) . corresponding to (11.31) . If le^u is the rate function for Z^'^ = (x — Xf) ■ e (i.e., 
lim^^oo^"^ lnPQ(Zf'^ > ct) = -Ie,u{c) for c> and any x), then Ie^u{cl{u, /)) = /'(O) holds 
for KPP reactions. From [25j we know that in two dimensions we have c*(u, /)/ \J De{u) = 
2V/'(0) +0(f (0)3/4) gj^g^u j/(Q) j^pp ^i^j^ (e, M)-independent error bound. 
This means that Ie,u{c) ~ c^/ADe{u) for c < a/ De{u), in the sense of 4:De{u)Ie,u{c)/c'^ — > 1 
as c/ ^JDf,{u) 0, uniformly in e,u. That is, effective diffusivity De{u) yields a good 
approximation of the rate function Ie,u{c) for c < ^/DJJi) in two dimensions. However, 
Theorem 15.11 shows that this is not true for general flows in more dimensions. At this point 
we do not know how to explain this difference. 

Motivated by a conjecture from [HIS] (see Corollary 1 1 . 3 1 b elow) . we are particularly inter- 
ested in the strong flow asymptotics of the unique/minimal speed cl{Au, /) for 

Tt + Au{x)-VT = AT + f{T), (1.10) 

with the flow profile u as above and flow amplitude A G M large. A natural question is 
which flow proflles are able to arbitrarily speed up fronts provided their amplitude is large 
enough (see [H [U [131 [13 (SSI [251 [23 [3Ql [32]). Having proved Theorem II. H this question in 
two dimensions becomes equivalent to the question of the asymptotics of De{Au). The latter 
is much simpler since De{Au) can be computed via (11.51) and (II. 6p with Au in place of u. 
In particular, it has been proved in [6l [121 [25] (see, e.g.. Proposition 1.2 in [25]) that in any 
dimension, limsup^^o^ i5e(^w) < oo when the equation 

M-V0e = M-e (1.11) 

has a solution 0e ^ H^iT'^), and lim^^oo Df.{Au) = oo when (11.111) has no such solution. We 
therefore obtain the following characterization. 

Corollary 1.2. Let u{x) be a 1-periodic incompressible mean-zero C^'^ flow on R^, let e G 
be a unit vector and f any combustion-type reaction, 
(i) If ([LTTD has a solution (pe e H^iT^), then 

limsupc*(ylu, /) < oo. (1.12) 

(a) If (II. lip has no solutions in H^iT"^), then 

lim cliAuJ) = oo. (1.13) 
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Remark. This result has been proved in [25] for KPP reactions but is new for general /. 

Of particular interest have recently been both percolating and cellular flows. Percolating 
flows possess streamlines (solutions of the ODE X' = u{X)) joining x-e = —oo and x-e = oo, 
a special case being shear flows u{x) = v{x2, ■ ■ ■ ,Xd)ei. Existence of such streamlines has 
obviously a strong effect on speed-up of fronts. Cellular flows, on the other hand, possess 
only closed streamlines, a prototypical example being Ucen{x) = V"'"(sin27rxi sin 271x2) whose 
streamlines are depicted in Figure [H Their effect on speed-up of fronts is therefore more 



Figure 1. Streamlines of the cellular flow Wceii- 

subtle, with diffusion across a thin boundary layer near the flow separatrices playing an 
important role. The interest in these flows stems from them being ubiquitous in nature. 
They appear as a result of instabilities in fluids such as Rayleigh-Benard instability in heat 
convection, Taylor vortices in a Couette flow between rotating cylinders, or heat expansion 
driven Landau-Darrieus instability. 

Corollary 11.21 shows speed-up of fronts in the sense of (11.131) for both percolating and 
cellular flows but known estimates on the effective diffusivity and Theorem 11.11 yield more 
precise asymptotics. For percolating flows 

< liminf ^^^^^ < hmsup^^^^^ < 00 (1.14) 

has been conjectured in and later proved for all / in |T7], and can also be recovered from 
Theorem 1.2 in [32] and Theorem 11.11 The asymptotic cKAuceii, f) ~ A^^'^ for the above 
cellular flow has also been conjectured in [D, [3] and obtained in [23] for KPP reactions, but 
the best result for general reactions has been A^^^ < cKAuceii, f) ^ A^^'^ for e = ei [T7] . 

Using Theorem 11.11 and the estimate D^^Au) ~ A^/^ from [2D], we now obtain the con- 
jectured asymptotic from [H |3] for general incompressible periodic cellular flows and all /. 
We also need to assume, as in [20], that the stream function of the flow has only non- 
degenerate critical points (otherwise the result is not true in general). That is, there is 
a periodic C"^'^ function : ^ M with only non-degenerate critical points such that 
u = V^H = {—Hx2i Hxi) and the complement of the the level set if = has only bounded 
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connected components {flow cells). Notice that this allows for almost arbitrary periodic 
geometry of the flow cells. We will call such u periodic non- degenerate cellular flows. 

Corollary 1.3. Consider a 1-periodic non-degenrate cellular flow u on M? , let e G he a 
unit vector, and f any combustion-type reaction. Then 

< limmf —- — < limsup --7- — < 00. 1-15 

As mentioned above, we also address the closely related question of quenching of reaction by 
cellular flows. Quenching occurs in the Cauchy problem for (11. ip with initial data T[0,x) = 
To(x) when \\T{t, ■)\\oo —>■ as t —>■ 00. This happens for ignition reactions when Tq is small 
in some sense so that T{t, ■) becomes uniformly smaller than the ignition temperature for 
some r > (and thus T solves (11. 3p for t > r by the maximum principle). But sometimes 
even large initial data can be quenched with the help of enhanced diffusion due to mixing by 
strong flows. 

It has been proved in [TT] that the flow Aucew quenches initial data supported in a strip 
of width L provided the amplitude A > L'^lnL. The authors also conjectured that the 
factor InL can be removed (heuristically, the minimal quenching amplitude should be such 
that cl{Au, f) ~ L, i.e., A ~ L^). We prove this conjecture for all periodic non-degenerate 
cellular flows on M^, which are also symmetric across the X2 axis (i.e., the stream function H 
is odd in Xi). We will call such periodic non-degenerate flows, which include Wceii, symmetric. 

Theorem 1.4. Consider a 1-periodic symmetric non-degenrate cellular flow u onM? . There 
are 7e > (independent of u) and Cu,e > such that if f is an ignition reaction with 
ignition temperature 9 > and ||/(s)/s||oo < 7e; then initial data Tq{x) G [0, 1] supported in 
[—L, L] X R are quenched whenever A > Cu.eL'^- 

Remarks. 1. We note that the bound 70 on / is in fact necessary. It is easy to show that 
if / is large enough, then even a single cell with Dirichlet boundary conditions can support 
the reaction for any A pT]. 

2. As in [11], scaling shows that ||/(s)/s||oo < 7e can be replaced by the requirement that 
the period of u be smaller than ||/(s)/s||^^7g 

The rest of the paper is organized as follows. In Section [2] we introduce our main technical 
tool. Theorem 12. H relating short and long time diffusivity of (11.31) in any dimension. In 
Section [3] we prove Theorem 13. H a relation between the speed and the width of a front in 
two dimensions, and then Theorem 11.11 In Section H] we prove Theorem 11.41 and in Section [5] 
we provide a counterexample to Theorem 11.11 in three and more dimensions. 

It is a pleasure to thank James Nolen for useful discussions. The author has been supported 
in part by the NSF grant DMS-0901363 and by an Alfred P. Sloan Research Fellowship. 

2. Short Time Diffusivity of Periodic Flows 

In this section we show that there is a close relation between short- and long-time diffusivity 
of the parabolic operator in (II. 3p . The results contained here are valid in any dimension. 
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Consider the stochastic process Xf starting at x G M'^ and satisfying the stochastic differ- 
ential equation 

dX^ = V2dBt- u{X^)dt, X^ = X, (2.1) 

where Bt is a normahzed Brownian motion on with Bq = (defined on a probabihty space 
{fl, Boo, IPq))- By Lemma 7.8 in [24], we have that if ip solves (11. 3p with ip{0, x) = ipo{x), then 

^{t,x)= [ My)kt{x,y)dy = En{MXn), (2-2) 

with kt{x,y) the fundamental solution for (11. 3p and Ej^ expectation with respect to u E Q. 
That is, kt{x, .) is the density for . 

We also have that if T solves (11.11) with T(0, x) = Tq{x) and ipo{x) = To{x) > 0, then by 
the comparison principle [26], for all t, x, 

< ip{t, x) < T{t, x) < e*ll^(')/'l'-V^(t, x). (2.3) 

We will use this relation to obtain short time upper estimates on the solution of (11.11) . 
We start the study of short time diffusivity for (11.31) with noting that uniformly in x G M'^, 

lim Eo (^^^^f^) = Deiu). (2.4) 

This is based on the fact that 

En{\{X^ -x)-ef) = [ \{y - x) ■ e\^kt{x,y) dy = [ \{Vty - x) ■ e\V^kt{x,Vty) dy 

and the following estimates on kt{x, y): 

f^^^ktix^Viy) ^ kl{0,y) in L'^{W'-) as t oo, uniformly in x; 

< kt{x, y) < Ct-'^/2e-l^-?'l'/^* for some w-dependent C > 0. 
The first is the standard homogenization limit (see, e.g., [HI [2T]) with 

k*ix V) = — g-(s/-a:)--D(«)-Hs/-^)/4t 

v/det (/}(«)) (47rt)^/2 

the fundamental solution of (11. 4p . the second is the Nash-Aronson estimate for mean-zero 
flows (see, e.g., [22]). The limit (12.41) now follows from 

The main result of this section is a lower bound on short time diffusivity for (11.30 : 

Theorem 2.1. There is C > such that for any r >1, any 1-periodic incompresible mean- 
zero Lipschitz flow u and any a > there are x G M'^ and t G [0, r] such that 

(|(Xf - x) ■ e| > aVrDe(M)) >l-Ca. (2.5) 
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Proof. We first note that if i7 = T"^ x f2 is equipped witli tlie product probability measure, 
tlien (12.41) gives 

lim [^^^^r^) = ^^-^^ 
wliere tlie expectation is witli respect to (x, u) G Cl. Tlien we liave 
Lemma 2.2. There is C > such that for any r > 1 and u as in Theorem \2.1\ 

¥.^{\{X:-x)-e\')>CrD,{u), (2.7) 

Proof. Let us first assume r = 1. For x G T"', let Xf = modi G T"' be the process 
corresponding to (11.31) on T'^. We note that Xf is uniformly distributed over T'^ as a random 
variable on f2. Indeed, if i? C T'^ and ipt){x) = Xb{x), then for each t >0, 

F^{X^eB)= I tlj{t,x)dx= I tlj{0,x)dx=\B\ 

because the evolution (11.31) preserves the total mass of if). 

We next let = Xf — Xf_^ and = ■ e. Then periodicity of u implies 

Y-X \A^t — 1 \A^t-l s^x 

t — ^1 - — ^1 - 

in law. Since the X^_-^ for alH > 1 are identically distributed as random variables on Cl, the 
same is true for the increment displacements Y^^ as well as the Z^. In particular, for each 
t > 1, 

E^i\Z^\')=E^i\Z^,\'). (2.8) 
We also have (X^ — x) ■ e = J2n=i N eN and so by fl2.6p . 

N 

{ZIZ"^) = ^ND^iu) + o{N). (2.9) 

n,m=l 

Since |Ej^ (Z,^'Z^)| < E^(|Zfp) is obvious from (12. 8p and the Schwarz inequality, we will 
obtain (12.71) for r = 1 if we can show the existence of w-independent M G N and 7 > such 
that 

|En(^:^i^)l<2-^l"^-"lE^(|Zrn. (2.10) 

whenever \m — n\ > M + 1. 

We denote ht{x, y) = J^jez^ ^t{x, y + j) the fundamental solution for (II. 3p on T'^. We then 
have 

ht{x,y)dy= / ht{x,y)dx = 1. (2.11) 



By Lemma 5.6 in there is M G N such that for all u as above, \\hM{-, •) — l||oo < |- The 
maximum principle gives ||/ij(.,.) — l||oo < ll^s(-)-) " l||oo for t > s and this together with 
{ht - 1) * {h, - 1) = ht+s - 1 (from (I^TTD ) implies for all m > M and 7 = (2M)-\ 

||/^„^(.,.)-l||oo<2-^("^+^). (2.12) 
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As a final prerequisite, we note that 



(2.13) 



This can be obtained by taking the initial datum iPq{x) = XT'* (2^) iii f!l-3p on M*^ and evaluating 
d 



dt 



I xijj dx = xAil) — xV-{uil))dx= I wp dx = / u{x) y ^ ^jj{t, x + j) dx = 



where we used integration by parts, the fact that Xljez^ = 1 ^.nd u being mean- zero. 

2' 2J 



Thus for each t > (recall that = [-i, ii*^^ 



= / xil){t^x)dx= / xkt{x,y) dxdy = / [x — y)kt{x,y) dxdy 

Jr^ jR^xTTd JlR'^xT'' 

because f^^ ykt{x, y) dx = y. Then by periodicity, 

^ / E + «)M. + , y) <^.<*. ^ - / E (i- - - « - 

In what follows we denote Xe = x-e for x G M'^. For m — n > M+l we have by fej^/cs = h+s, 



[X'fi ^n— l)e (^m ^m— 1 )e^n— 1 (^5 ^n— 1 ) ^1 (^n— 1 5 ^n) 



hjji_\_^i{xYi'i X'fYi—\^k\{xYn—\^ x^n ~t~ ^) dxdX'^—idxYidXjYi^—idx^^ 



Here we have used the fact that 

g, + / + g) 

^ ^ (-^m ~l~ ^ -^m— 1 )e^l (•^m— 1 j -^m ~l~ ^ ^ ^m— 1— n (-^n ; "f'm— 1 ~l~ •?) 
^ ^ (-^m ~l~ ^ Xm,~l)ehm~l—n{Xni Xm—l)k\{Xm—l^ Xm ~\~ 



igZd 
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(due to periodicity of u) and similarly 



^ ^ ("^n ~l~ J ■Xn—l)ehn~l iS^i ■^n—l) ^1 (-^Ji— 1 5 -^n ~l~ J ) • 



The integral with respect to x can now be eliminated along with hn-i{x,Xn-i) because 
hn-i{x,Xn-i) dx = 1. Notice that fl2.13p gives 



/ y^X^m + I - 



■Xm—l)ek\{Xrn—l^ •^m ~\~ dXm—ldXm 0, 



and so (12.121) and Schwarz inequality imply 



2'm— 1 )e^l (-^n— 1 ; ~l~ j) 



[^m— 1— ?i('''n) '^m—l) l]^l('^m— 1; "^m ~l~ dXn—\dXndx^—\dXu 



< 



1/2 



1/2 



1/2 



Now (12.81) yields (12.101) . finishing the proof for r = 1. The general case is identical, this time 
with being X^^ — X1^^_-^-^^, the same 7, M and (7, and 2N replaced by 2Nt in ( 12. 9p (one 
actually gets C* — > 2 as r — > 00). □ 

We will now prove (12.51) with C = lOC*^^/^. Assume, towards contradiction, that for some 
u there are r > 1 and a > such that for any x G M*^ and any t G [0, r]. 



Pf, (\{X^ -x)-e\< a^/rDeiu)^ > Ca. 



(2.14) 
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We first claim that for eacli x G M*^, 



VtG [0,r] 



\{X^-x).e\< 



c 



> 



(2.15) 



Indeed, if this is not true, let x be such that the probability in fl2.15p is less than i. This 
means that there is a subset fl' <Z fl of measure more than | such that if tj^uj) > is the 
first time the (almost surely continuous in t) path = Xf{uj) hits the set 



H, = {yER''\\iy~x)-e\= 2ja./^Dj^}, 

then tj{uj) < T for each uj E Q' and j = 0, . . . , [^J- Now the strong Markov property of 
the process Xf , the fact that t — tj G [0, r], and fl2.14p imply that for j = 0, . . . , [^J the 
following conditional probability satisfies 

Fn {\iX^ -x)-e\e ((2j - l)a^/^D^, (2j + l)«VrDe(M)) | i?^^,) > Caxnicu), 

with Btj the cr-algebra corresponding to the stopping time tj. Thus 



Fn {\iX^ -x)-e\e ((2j - l)a^rD,{u), (2j + l)a^TD,{u) 



> 



Ca 



for j = 0, . . . , [^J- Since ^(Lc^J + 1) > this is a contradiction, thus proving (12.151) . 
Now (I2.15P and the almost sure continuity of in t show for each x G M'^ and each j > 1, 



Fc 



\{X:-x)-e\> 



C 



< 



That, however, means 



{XI 



X 



2 



96rDe(M) 



< CtD,{u) 



contradicting (12. 7p . This finishes the proof of Theorem 12.11 



□ 



3. Pulsating Front Speed in 2D 

We will now prove Theorem 11.11 The upper bound in (11.71) is immediate from the results 
in [25j. Indeed, the same bound has been proved there for KPP reactions with C2{f) = 
C a/ /'(0)(1 + a/ /'(O)). One thus only needs to apply this result to some KPP / such that 
/ < / and ||/(s)/s||oo = /'(O), and use the fact that c*(m, /) < c*{u, f) due to / < /. 

We will now prove the lower bound in (II. 7p using Theorem 12.11 coupled with the following 
bound on the width of the front in terms of its speed. We note that we can assume Tt > 0. 
This has been proved in [2] for all ignition reactions (and also for positive reactions with 
/'(O) > 0), and the lower bound in (11.71) for any ignition reaction / < / proves the lower 
bound for / because again c*(m, /) < c*(-u, /). 
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Theorem 3.1. There is Cq > such that for any 1-periodic incompressible mean-zero C^'^ 
flow u on M^, any unit vector e G M^, and any combustion-type reaction f the following holds. 
If f ^ ^XlcQ for some m > and < C < ^ < 1, £^ G (0, (^ — C)/2), and T(t, x) is a pulsating 
front for (11. ip with speed c > and Tt > 0, then there is z & M. such that any connected set 
B with 

B (Z{xe^^\x-e>z + ct + cCo{m~^ + e'^) + 2 and T{t, x) > ( + e} = B+ 

or 

BC{xeR'^\x-e<z + ct- cCo{m'^ + e'^) - 2 and T{t, x) < ^ - e} = B' 
satisfies diam(i?) < jq. 

Remark. The above form of this resuh will be sufficient for our purposes. Its proof in 
fact shows that the set of x such that T{t,x) > C + resp. T{t,x) < ^ — e covers at most 
1% of any unit square lying in the halfplane x ■ e > z + ct + cCo{m~^ + + 2 resp. 
X ■ e < z -\- ct — cCo{m~^ + — 2 (and this bound decreases as Cq increases). That is, for 
any t, outside of a strip of width 2cCo(m^^ + e~'^) + 4, values of T{t, ■) are mostly outside of 
[( + e,C, — £]■ This yields a bound on the width of the front in terms of its speed. 

Proof. Let T(t,x) = U{x ■ e — ct,x) with U satisfying (11.21) . so that 

-cUs + u-VxU + u- eUs = A.,U + f/,, + 2e ■ V,f/, + f{U). (3.1) 

Integrating this over F = M x Fq = M x [0, 1]^ and using 1-periodicity of U in x, (II. 2p . and u 
being incompressible and mean zero, we get 

j f{U{s,x))dsdx = c. (3.2) 

Similarly, multiplying (13. ip by U and integrating over F yields 

- + / {V xU + eUsl'^ dsdx = / f{U{s,x))U{s,x)dsdx. 
2 Jr Jr 

The right hand side is bounded above by c thanks to (13. 2p and so 

J f{T{t,x))dtdx = 1, (3.3) 

/ \VxT{t,x)\^ dtdx < -. (3.4) 
Jr 2 

Let T(t) = T(t, x) dx G [0, 1] so that Tf > because Tt > 0. Assume that T(t) G [C + 

|,^ + |] for some t. Then the Poincare inequality \\T(t, ■) — T'(t)||2,2(ro) < Ci\\'VxT{t, ■)||L2(ro) 

for some Ci > 2 gives either ||r(t, ■) - f(t)||i2(ro) < | or \\VxT{t, ■)\\L2(^ro) > 4§7- In 

the first case the set of x G Fq with \T(t,x) — T{t)\ > | has measure less than ^ and so 

j^J{T{t,x))dx > f . Thus (EJl) and (I33D show that the interval of all t with f{t) G 

[C + |, ^ + f ] has length at most ^ + 

Assume now T(t) < C + | and that there is a connected set B <Z {x E Tq \ T(t, x) > C + ^} 
with diam(i?) > ^ (such t form an interval due to Tj > 0). Projection of B on one of the 
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axes (let us say Xi) is then an interval / of length more than ^. If for each xi G / there is 
X2 G [0, 1] such that r(t,xi,X2) < C + T' ^^^^ \\^xT{t, ■)llL2(ro) ^ Jq- If the other hand 
T(t,y,X2) > C + ^ for some ?/ G / and all X2 G [0, 1], then there is J C [0, 1] of measure at 
least I such that either T(t, Xi, X2) > ^ ^^r all (xi, X2) G [0, 1] x J or T(t, Xi, X2) < y for each 
X2 G J and some X2-dependent xi G [0, 1]. In the first case \\T{t, ■) —T{t)\\L2(j'^^ > |, and the 
Poincare inequality gives || Va;T(t, ■)||i2(r(,) > In the second case || Vxr(t, •)||x,2(rg) > j^. 
In either case we have ||Va;T(t, ■)||L2(ro) ^ so the interval of the t above has length at 

, 32C2 

most — 

The same is true for t such that T(t) > ^ — | and there is a connected set S C {x G 

To I T(t, x) < — £:} with diam(i?) > ^. Thus there is an interval [a, b] with b — a< ^ + 
such that for t < a connected subsets of {x G Fq | T(t, x) > ( + e} have diameter at most 
and the same is true for t > 6 and connected subsets of {x G Fq | T{t, x) < ^ — e}. 

Finally, let z = — and Cq = 36C^, and assume that for some t a set B C Bf contains 
a point X. Then with [xj the integer part of x, 

T(t,x) = U{x-e-ct,x) = U{x-e-ct,x- [xJ) =T(t-^,x- [xJ). 
We have using x E B, 

t _ \lLl < ^ _ f _ Co(m-i + e-') --< ^ - 36Cf{m-' + 5-^) < a 
c c c c 2 

and so diam(i?) < ^ (in fact, to handle B not lying entirely inside a square with integer 
corners, we need to consider Fq = [— ^, 1 + which only changes Cq by a fixed factor). A 
similar argument takes care of sets B C B^ , thus finishing the proof. 

Notice that if B is not required to be connected in the proof, then one still obtains the 
bound jg on the size of its projection on the axes. This proves the remark. □ 

Proof of Theorem As mentioned above, we only need to prove the lower bound in (11. 7p 
and only assuming Tt > 0. For ( G (0, 1) let ^ = 1— and e = , and pick M^, > 
so that with C from Theorem |2JJ 



-M, 



^-e 1-iC + e) 



< 1 - Ca^. (3.5) 



This is possible because if C' > C + ^ and 1 — ^'^ ^ < ^ — e (e.g., (' = ^-^), then 

e-e i-(C + ^) i-il^^i-C 2(i + 2C-C^) • 

Let us choose ( so that / is strictly positive on [(,(,] and denote m = min^g^^ ^] /(s) > 0. It 
is sufficient to assume that 

m 

m < MfC and f{s) < —s for s G [0, 1]. (3.6) 
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Indeed, otherwise consider /(s) = min{/(s), ^s, M(^s} (which satisfies (13.61) with rh = 

minsg[^^g] /(s) = min{m, M^^} in place of m) instead of / and then use cl{u, f) < cl{u,f). 
Note also that ([L9D for all m^(/) G (0, M^C) proves ^M) for all m^(/) > (with a different 
C^). So let us assume (13. 6p . 

Let T be a pulsating front for (11.11) in direction e and with speed c*(-u, /), assume z = 
in Theorem 13. II (otherwise shift T in time by z/cl{u, /)), and let r = Mc_C,/m > 1. Theorem 
12.11 and (13. 5p show that there is x G and t G [0, r] such that either 

C+£ 



P.c 



((X,^ - x) ■ e < -acv^^^D^^) > ^ (3.7) 



or 



(e - ^)e 



(^(Xf - x) ■ e > ac V^^e(«) j > ;l ^ • (3.8) 

We now claim that 

a^^/rD.iu) - ^ < c:(m, f)t + 2cl{u, f)Co{m-' + e~^) + 8. (3.9) 
If so, then t < t = M(^(/m gives 



Ci^/mD,{u) -C2m 
63 + 64m 

with some (^-dependent positive constants. Since De{u) > 1 by (II. 6p . this then gives 
C5^/mDe{u) < cl{u, /) for some C5 > and all small enough m > 0, thus proving (ll.7p and 
m for all m^(/) > 0. 

It remains to prove (13.91) . and it is sufficient to consider ac^\jTDf,{u) > -j^. Assume first 
(13.71) . Due to spatial periodicity of u we can assume for x,t from (13.71) . 

\x-e- cliu, f)t - cliu, f)Co{m-^ + e"') - 4| < 1. (3.10) 

We can also assume 

T(t,x)<C + £ (3.11) 

at the expense of changing (13. 7p to 



¥n ((Xf -x)-e<^- a^^/rD^iu)) > |— ^. (3.12) 

This is because Theorem 13.11 shows that the largest connected set B of points satysfying 
(I3.10p but not (13. lip and containing x (if it is non-empty) has diam(i?) < (recall that 
z = 0). Hence almost sure continuity of in t shows for any P < —j^, 

Fn {{X^ - x) ■ e < /3) < sup ((Xf - x) ■ e < 13) < sup Pf, ((X,^' - y) . e < (3 + ^) . 

yedB yedB 

se[o,t] se[o,t] 

Let r be the union of all connected components of the set F' = {y G | T(0, y) < ^ — e} 
lying entirely in {y E M."^ \ y ■ e < —c*{u, f)Co{m~^ + e~'^) — 2}. Then by Theorem 13.11 each 
connected component of F has diameter at most ^ and 

F D r n G I y ■ e < -c:(m, /)Co(m-i + e'^) - 3}. (3.13) 
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Let ip solve (11 .Sp on the domain \ F with '?/'(0, y) = T(0, y) for ?/ G \ F and ip{s, y) = 
T{s, y) {>^-e because Tt > 0) for y e dT. Then similarly to ([22D and flO) . 

ij{t, x) = En (T(0, Xf)xa>t + Tit - a, X:)x.<t) > En (T(0, Xt)Xa>t + - e)xa<t) , 

(3.14) 

with a = <7{uj) = infxf ((^)Gr s, and 

< ipit, x) < Tit, x) < e*'"/'^V(^, x). (3.15) 
If (E3D is violated, then fIXTU]) and flXT^ show that 

(Xf ■ e < -<(«, f)Co{m-' + 5-2) - 3) > 

so that fl3:T3|) and fIXTD yield ^(t,x) > |^(^ - e) = ( + e. But this contradicts flXTTD and 
(I3.15p . so (13. 9p is valid and we are done in the case when (13.70 holds. 

Let us now assume (13.80 . Here one uses a similar argument with (13.100 . (13. lip , and (I3.12p 
replaced by 

\x-e + cl{u, f)Co{mr^ + e~'^) + 4| < 1, (3.16) 
T(0,x)>^-e, (3.17) 

¥n[{X: - x) ■ e> a^./^D:(uj - ^) > ^ _ ; (3-18) 

and F the union of all connected components of the set F' = {y G | Tit, y) > C + ^} lying 
entirely in G | ?/ ■ e > c*(m, f)t + cHu, f)Coim~^ + e^^) + 2}. Again each connected 
component of F has diameter at most ^ and 

T n {y eR^ \ y ■ e > cliu, f)t + cHu, f)Coim-' + e-^) + 3}. (3.19) 

If dMD is violated, then fl3A6|) and fl318D show 

1 — (f — rV^*™/'^ 

{X: ■ e > cliu, f)t + cliu, f)Coim-' + e-') + 3) > | , 

and so (13.191) and 

ijit, x) = En (r(0, Xnx.>t + Tit - a, X:)xa<t) < En (T(0, X^X.yt + (C + e)Xa 
(with ip defined as above) yield 



i-(C + ^) ' V i-(C + 

This again contradicts (13.151) . (13.171) . and Tt > 0, so (13.91) is also valid when (13.81) holds. □ 
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4. Quenching by Cellular Flows 

Proof of Theorem \1.4\ It is obviously equivalent to consider the problem on M x T with Tq 
supported in [—L, L] x T, which is what we will do. It has been proved in [20J that periodic 
non-degenrate cellular flows in two dimensions have Df,(Au) ~ A^^'^ as A — > oo (the special 
case of Mceii has been treated earlier in [12]). It is therefore enough to prove that a 1-periodic 
symmetric non-degenerate cellular flow u quenches solutions of (11. ip when 

\\f{s)/s\\^<^g (4.1) 

and To is supported in [—bg ^J De{u), be Deiu)] x T, with e = ei = (1,0) and some be > 0. 
Having Theorem 12.11 at hand, the proof is similar to that of Lemma 4.2 in [30] . 

Let us prove the last claim. We only need to consider large enough De{u) because solutions 
supported in [—be,be] are quenched for any u provided (14. ip holds and 6^ > is small 
enough [2S1 Theorem l.l(ii)]. By the comparison principle, it is sufficient to consider initial 
data Tq{x) = X[-L,L]ixi) with L = [■^^y De{u) — |J and C from Theorem 12.11 (then we 
can take be = 6^(1280)-^ mm{l,2b'g/ 3} and have be^/ojuj < max{6'g,L}). Let ■0 be the 
solution of (II. 3p with initial datum i/jq = Tq. We first claim that there is a continuous curve 
h : [0, 1] ^ [0, 1] X T such that (/i(0))i = and (/i(l))i = 1 , and for all s e [0, 1] and t > 1, 

^{t,h{s))<^-. (4.2) 

To this end we let cf) be the solution of (11.31) with initial condition {/)(0, x) = Xi~k-2,k]{xi) 
where K = [8L6'~^'|. Theorem 12.11 with 

K + 2 8L + Ae 9 
a = —^== < == < — 

,/Dju) - e^/Dj^ - sc 

shows that there are r < 1 and y G [— 1, 0] x T such that 

0(r, y) = Fn{{Xy)i E [-K -2,K])<^-. 

The maximum principle for (11.30 implies that the connected component of the set 

{{t, x) G [0, r] X M X T I 0(t, x) < |} 

containing (r, y) must intersect 

{x G M X T I 0(0, x) < |} = (M \ [-K -2,K])xT. 

Since by symmetry (f){t, Xi, X2) = (pit, — 2 — xi, X2) for xi > 0, this means that there is a curve 
h{s) joining {0} x T and {K} x T such that for each s there is < r with 

0(r„ his)) = ¥n ((X^(^))i G [-K - 2, K]) < ^. (4.3) 
Lemma [2.1( iii) in |30j shows for any M G N, t > 0, and any 1-periodic symmetric fiow, 

\yi\ 

M 



Fn{\Xy\<M) < 
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Applying this with M = L, y = X^J'*'' and t — Tg (for t > 1 > r^) in place of t, we obtain 
using (14. 3 P for any t > 1 and s, 



9 





'K' 


-) 

8/ 


T 



<.+ (!-.).<- 



which is (14 ■2p (after reparametrization of h and restriction to s G [0, 1]). 

Symmetry of u and tpo implies that (14 ■2p holds for h{s) extended to s G [—1, 1] by h{—s) = 
{—{h{s))i, {h{s))2). Finally, (14.20 applies to h{s) extended periodically (with period 2) onto 
M. The last claim holds because ip{t,x) > ilj{t,x + (2,0)) when xi > —1 (and ilj{t,x) > 
ip(t,x — (2,0)) when xi < 1). This in turn follows because ipit^x) — ipitjX + (2,0)) solves 
(ll.3p with initial datum that is symmetric across xi = —1 and non-negative on [—1, oo) x T, 
and hence stays such by the symmetry of u across Xi = —1 (the latter is due to the symmetry 
of u across Xi = and periodicity). 

This means that \\ipit + 1, ■)||oo < Wvi't, OIloo + f where r] is the solution of (ll.Sp on the 
torus [—1, 1] X T (with —1 and 1 identified) with ri[0,x) = 1 and r]{t,h{s)) = for t > 
and s G [—1, 1]. Then Lemma 2.3 in [30] shows that there is a universal constant 6 > such 
that ||^(t + l,-)||oo < 2e-^* + f. We let = ^ In f + 1 so that \\ipiTe, ■)\\oo < f- If now 
7e = Tg^ In 2, then 

||T(re,-)||oo<e^''"1l^(re,-)||oo<^. 
The maximum principle then implies that for t > tq the function T solves (II. 3p with 
\\T(t, ■)||oo < 0, and quenching follows. □ 

5. A Counterexample in 3D 



In this last section we show that our main result, Theorem ll.il does not hold in three and 
more dimensions. The counterexample we provide also shows a breaking of symmetry in more 
dimensions. Specifically, despite the fact that Df.{u) = D_(.{u) for any e and m, it is not true 
that c*(u, f)c*_^{u, f)~^ is bounded by w-independent positive constants when d> 3. 

Theorem 5.1. Assume d > 3, e = ci, and f is any KPP reaction. Then there is no Ci > 
such that for any 1-periodic incompressible mean-zero C^'^ flow u on M*^ we have cl{u, f) > 
Ci^y De{u) , and there is no C2 < oo such that for any such u we have c^iu, f) < C2a/^ 

Proof. We will use the following two results from [32]: 

cliAuJ) J^^{u-e)w^dx 
lim = sup — — 7-^ . 

2< " 



||V«;||2</'(0)|hll2 



JDe{Au) La{u-e)wdx 
lim — = sup , 

where 

Z= {w^H^iJ'^) I M-Vw = 0}. 
We let e = ei and we will assume d = 3 since the general case is analogous. Let x ^ be 
a smooth characteristic function of the unit disc in and for small i? > let XR^P^^^'^'i) = 
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j,k€i? x(^V^5 ^^)- ^"^^^ consider the mean- zero periodic shear flows Aur{x) — 
(A'Uk(x2, X3), 0, 0) with A G M and ur{x2,x^) = Xr{.^2)X^) — 1- In this case the elements of 
X are precisely the w G H^(T^) which are independent of Xi and the above formulae become 

liin ^*e{^^RJ) ^ /t2 UrW^ dX2dX3 

iiv«.iii</'(o)ihiii 

hm <i^!^= sup U--nWdx2dx,^ ^^^^^ 
\A\ ^(zx> \\w\\2 



l|V«>|ji</'(0)||«; 

^-±~ |A| ^6j/ ||Vti;||2 

where T' = H^{T'^) and e = ei. 

For (x2,X3) G [-|, = andr = a/x| + x| let iyij(x2,X3) = max{0, minjlog 27^, log ^}}. 



Taking w = ty^j we see that RHS of (15.31) > ylog for small R. On the other hand, —ur < 1 

and so RHS of ( 15. 2p < 1. This proves the first claim. 

Now let wr be the maximizer of the RHS of (15.31) . which exists by [32], normalized by 
||Vwr||2 = 1 and fr^2U!iidx2dxs = 0. The Poincare inequality then gives ||wi?||2 < C for 
some C > 1. Let C = maxjC, (/'(O))-^}, choose Kr e R so that \\wr + Kr\\1 = C, and 
define wr = wr + Kr. Then 1 = ||Vwi?||2 = Ilw^ijlli/C* < /'(0)||w/?||2 (so wr enters in the 
RHS of (15. ip ) and wr also maximizes the RHS of (15. 3p . Schwarz inequality, IlifijUl = C, 
Ur = Xr ~ ^1 &nd the previous paragraph then imply 



(5.4) 



/ XRWRdx2dxz + \/C > / URWRdx2dx3 = RHS of (EJD > J^ogj^ 

for small R. But then taking w = wr in (15. ip and using Schwarz inequality, Jj2 Xr dx2dx3 = 1, 
and (15.41) gives 

. i jj2XRWRdX2dX3j 

RHS of (EH) > i / XRwldx2dx3 - 1 > ^ ^ ^ - 1 > (RHS of (ED)^ 

Jt2 ^ Jj2XRdX2dX3 

for small R. This together with (15. 4p proves the second claim and we are done. □ 
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